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SUMMARY 


An analytic method for the design of automatic controls is 
developed that starts from certain arbitrary criterions on the 
behavior of the controlled system and gives those physically 
realizable equations that the control system can follow in order 
to realize this behavior. The criterions used are developed in 
the form of certain time integrals. : 

General results are shown for systems of second order and of 
any number of degrees of freedom. Detailed examples for 
several cases in the control of a turbojet engine are presented. 
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INTRODUCTION 


In the past several years, there has been increasing deyelop- 
ment and interest in automatic control; in the fields of gun 
direction, guided-missile control, and control of gas-turbine 
engines, for instance, where very refined and accurate 
controlled behavior is required, need still exists for further 
development of the methods of controls analysis and design. 

Recent developments in this field have been coneerned 
chiefly with the problem of control analysis both in the 
realm of linear systems (reference 1) and in the realm of 
nonlinear systems (reference 2). These analytical works 
answer questions concerning the behavior of a given system 
or the effect on its behavior of changing certain constants 
in the system. 

Another problem of equal and, in some cases, greater 
importance is that of control synthesis. Work on this prob- 
lem seeks to find the method of design when certain cri- 
terions concerning the behavior of a controlled system are 
given. The method decided upon should give all aspects 
of the system; for instance, where the system should be 
linear, what the general configuration should he, and what 
the precise values of all the constants should be. 

This synthesis problem has hardly been broached in liter- 
ature. The use of analysis as a design procedure offers a 
partial solution to this problem in that the analysis of a 
large number of cases may reveal, coincidently, one that 
satisfies the desired criterions of controlled behavior. Such 
a method is, at best, long and tedious and almost always 
would result in compromises because the systems chosen to 
be analyzed would probably be such that they could never 
satisfy all the desired criterions. 

A method for designing a linear system to satisfy certain 
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special criterions when operating on a random input is 
developed in reference 3. This method is applicable as an 
addition to a control system whenever random external 
disturbances are involved. A partial solution to the syn- 
thesis problem is developed in reference 4 in satisfying the 
criterions of noninteraction for systems with many degrees 
of freedom. 

An analysis made at the NACA Lewis laboratory during 
the latter part of 1950 and presented herein develops a 
rational method of control synthesis that starts from any 
arbitrary but physically realizable criterions and results in 
the equations for the best system that satisfies these cri- 
terions. As is shown, the nature of the eriterions, in general, 
requires minimizing certain time integrals by using the 
calculus of variations and the methods developed are an 
application of the calculus of variations to the problem of 
control synthesis. 

A careful scrutiny is first made of the whole problem, fol- 
lowed by a development. of general results. These general 
results are then applied in examples to the design of turbojet- 
engine control systems. In general, the methods used vary 
according to the order of the differential equation describing 
the plant, the process, or the system being controlled and 
according to the number of degrees of freedom or independent _ 
variables. Detailed analyses are presented for application 
to a system of first order and of one degree of freedom. 
General results for systems of second order and any number 
of degrees of freedom are developed in the appendix. l 


SYMBOLS 


The following symbols are used in this report: 


a,b,c gas-turbine-engine characteristic constants 

C constant 

E function of A 

F,H functions of y and ġ 

Fa) indicates Ffy), y] 

F,(y, ý) indicates partial differentiation with respect 
to y 

ST So arbitrary function : 

Fide, - weighting functions used for gas-turbine con- 
trol criterions 

G, Gu G functions of y 

g. function used for gas-turbine-engine surge 


criterion 


upersedes NACA TN 2378, “Automatic Control Systems Satisfying Certain General Criterions on Transient Behavior” by Aaron 8. Boksenbom and Richard Hood, 1951. 
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h function used for gas-turbine-engine blow- 
out criterion 

K controller gain 

L, Li, E temperature limits for gas-turbine engine 

N, P, T, w; deviations of gas-turbine-engine speed, com- 
pressor-discharge pressure, characteristic 
temperature, and fuel flow, respectively, 
from values at some common equilibrium 
condition 

Ne actual engine. iccse 

n exponent 

t time 

i time at end of transient 

ôl: variation in time at end of transient a 

w, Y, 2, 

bw, y, 82 independent variables, functions of time 

€ small number 

A, An, Ag, Ag arbitrary constants 

G gas-turbine-engine time constants for re- 
sponse to temperature 

T transient time constant of controlled system 

Subscripts: 

(g initial condition of variable 

f final value of variable 

x setting or desired value of variable 

Superscripts: 

Ed 


indicates case different from optimum 


The dot indicates differentiation with respect to time. 
The prime indicates differentiation with respect to the argu- 


ment shown. 


ANALYSIS, 


SURVEY OF PROBLEM 


Control problem.—An important aspect of the control- 
synthesis problem is a clear definition of the criterions of 


desired controlled behavior. 


If a variable y is to be con- 


trolled, a reasonable criterion is that the time integral of 
some function of y is to be a minimum or a constant; that is, 


or 


$ 
í ! #(y) dt=constant or minimum (1) 
a 


t 
f ' (y— y} dt=constant or minimum (2) 
«70 


Equation (2), for instance, weights the error in y as the square 


and according to the time duration of that error. 


Another 


type of criterion may be that which requires a certain time, 
duration to be a minimum or a constant; that is, 


t 
[ ' dt=constant or minimum (3) 


F The use of a single criterion, such as equation (1), will 
‘usually yield f(y)=constant. This result is reasonable be- 
cause f(y) can usually be made identically a constant if no 
additional criterions are imposed on other variables in the 





system. Certain limiting conditions usually exist on other 
variables in the system and these conditions must be included 
in the original criterions. 

Thus, a possible criterion could be written 


t 
f ; (y—y,)? dt=minimum 
0 


a 


É 
T Fa) dt=constant 
0 


If, for instance, y=engine speed and z=characteristic tem- 
perature of a gas-turbine engine, the criterion of equation 
(4) states that it is desired to design a control system such 
that, for a particular value of a temperature integral, the in- 
tegral of the speed-error squared is a minimum. This eri- 
terion may be used if, for instance, it is known that an over- 
temperature condition ean be tolerated for a certain period 
of time and it is desired to keep the average speed error at 
a minimum during a transient. 

The general theory will show that as many criterions as 
desired of the type shown in equations (1) to (4) can be in- 
cluded together and a control system can be derived that 
automatically satisfies all these criterions simultancously. 

Another aspect of the control criterions is the end con- 
ditions of the integrals of equations (1) to (4). The time in- 
terval for which these integrals are to be a minimum or a 
constant must be chosen. A reasonable time interval is any 
duration during which essential external disturbances are 
constant and during which the system to be controlled moves 
from one essential levelof operation to another. The essen- 
tial external disturbances are those that cannot be imme- 
diately corrected by the control system. If an essential 
external disturbance were allowed within the time interval 
of the criterions, no physically realizable system could be 
expected to anticipate this disturbance so as to behave prop- 
erly before this disturbance occurs. An essential level of 
operation is any specific condition of only those variables 
that must be continuous. It will be shown that the essential 
level of operation appears as a natural boundary condition 
for the type of criterion used. In the case of a turbojet 
engine the transient behavior of which can be described by 
a first-order differential equation, the engine speed deter- 
mines the level of operation. If a lag exists in the fuel 
system or between temperature and engine speed, then both 
engine speed and acceleration are required to deseribe the 
essential operating level of the engine. 

Analytic problem.—The control system resulting from any 
design method must be physically realizable. There are 
two aspects to this problem. First, it is possible to set down 
eriterions that are not realizable with any system or are 
incompatible with each other. If such criterions are used, 
the unrealizability will appear either as a requirement on 
the control to look ahead into the future or as an inability 
to satisfy the boundary conditions of some differential equa- 
tion. In most cases, a clear understanding of the criterions 
used and of the system to be controlled will indicate incom- 
patibilities. of this sort. 
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The seeond aspect of physical realizability is purely mathe- 
matical. It is desired to derive a description (a differential 
equation) of the control or the controlled system that satis- 
fies the criterions of control and all the necessary boundary 
conditions that arise in the derivation of this equation. 
Although the mathematical solution of the problem may be 
any derivative or integral of this differential equation, the 
physical solution of the problem requires the differential 
equation that itself satisfies the boundary conditions and for 


which no undetermined constants of integration exist. Thus, 
such forms as 
y=C£ 
and (5) 
y=Cz 


are not necessarily interchangeable as descriptions of some 
part of a controlled system because the forms differ by an 
undetermined constant of integration. For stable linear 
systems, the effect of this constant becomes vanishingly 
small; for the general nonlinear systems presented herein, 
however, this constant must be considered. 

Stability problem.—The requirement of stability is a 
special criterion that does not enter into the main body of 
the methods of this report. It may enter in the final steps of 
the method where the final differential equation describing 
the controlled system may be the integral of a higher-order 
differential equation that satisfies the necessary boundary 
condition for stability. In addition, it is always necessary 
to add to the controlled system a stability device that does 
not affect the behavior of the system as far as satisfying the 
other criterions is concerned. This device can be described 
as follows: 


y=0 when y=y; (6) 
or, for a second-order system, 
y=0 when y=y, and y=0 (7) 
GENERAL THEORY AND RESULTS 


It has been shown that the criterions for control can be 
developed in the following forms: 


fy 
f f(y) di 


{ “y—y)?dt 
l À (8) 
ili f(a)dt 


foa 


and so forth. If, for such a list of criterions,-one of the 
integrals is to be a minimum under the condition that the 
other integrals are to be constant, it is sufficient (reference 5) 
to make 











t ty 
("jodin f uydi 


é t 
def Hadit | ' di=minimum (9) 
(03 0 7 
or 


li LAMA y—ys)? + dofo(2) +g] dé=minimum (10) 


The Ns are arbitrary constants that enter into the control 
system as the adjustable parameters and are precisely de- 
termined by the choice of values that the constant integrals 
are to have. 

The technique of the \ multipliers is widely used for prob- 
lems of this type where one condition is to be a minimum 
under other restrictive conditions. Indeed, the conditions 
need not be in integral form and any functional or differential 
relation among variables can be handled in a similar manner 
(reference 6). 

Equation (10) can be made very general when all possible 
restrictive conditions are included. In the final equations, 
which are derived later, if any one criterion is not to be used, 
then the corresponding A—>0. If any of the criterions is to 


be zero (as for the case of a variable having an absolute limit), 


then the corresponding \> œ. 
If, for this development, the system to be controlled is of 


first order and of one degree of freedom (has one independent __ 


variable), then the variables y and z are such that z=2(y,7). 
Equation (10) can then be written, in general, as 


é 
f Fly ý) di=minimum (11) 
Jo 


In equation (11), F is a continuous function of y and y, and 
y is a continuous function of time. ; 
The calculus of variations (reference 5) is used to determine 
y as a function of time such that the integral of equation (11) 
is a minimum; that is, if the solid curve of figure 1 makes 
the integral a minimum, any other curve (such as the dashed 
one) will make the integral equal to or greater than the first 


y YFKE) 








0 t t, 


t 
1 
FIGUBE 1.—IDustration of curves for minimization of Í Fiy, p dt. 
o 
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integral. If the two curves are very close, the condition 
that the two integrals are equal will make the integral of 
equation (11) stationary (maximum, minimum, or inflection 
point). 

The problem ‘is similar to that of finding the maximum or 
minimum point on a curve by setting the first derivative 
equal to zero. Whether a maximum or minimum point 
exists is decided by the second derivative at that point. In 
the variational problem, proving a true minimum involves 
taking the “second variation.” As in the problem of finding 
a minimum on a curve, the second variation proves, at most, 
a local minimum and not an absolute minimum. For the 
specific examples discussed herein, an absolute minimum is 
proved by another method. In many cases, the physical 
meaning of the equations will indicate the existence of a 
unique minimum that obviates going further than the 
“first variation.” 

On figure 1, the curve for minimizing the integral of equa- 
tion (11) is the solid curve having the value y(t) at any time. 
Any other curve differing from it by a small amount is shown 
as the dashed curve having the value y(é)+e y(t) at any 
time, where e is a small number and 6y(é) is an arbitrary con- 
tinuous function of time. The condition for the integral of 
equation (11) to be stationary is 

d fp +edey k : 
Ted Fly+edy, y+ dy) di=0 (12) 
at ` 
e=0 


where £ [dy (t)] = 6y (£) (reference 5). The time duration for the 


integral of equation (12) is such as to start at some definite 
time (t=0) but not to end at a definite time but rather 
along some curve y=f(é) (fig. 1), in order to allow the proper 
boundary conditions of moving from one essential level of 
operation to another, as previously discussed. 

Performing the operation indicated by equation (12) leads 
to 


t é 
. f ! Fy dy dt+ (Fs di+-F(t,) èti =0 (13) 

fo JO 
Integrating the second term of equation (13) by parts gives 


fi d n : 
Í |F, -4 œ) | èy arian] FG) = (14) 


Because ôy is an arbitrary function, the integral and the 


boundary-condition terms must vanish separately. The 
geometry at the end condition ¢=t, fig. 1) gives 
ôy =[F E) yE) èt . (15) 


as e>0. The two conditions that follow from equation (14) 
then become 


ti d 4 
| Í |r, — Fo) | èv dt=0 (16) 
and 


at [F-E 6) E — ýt) } —F; (0) dy(0)=0 (17) 








The time interval during which the criterion of equation 
(11) is to hold is considered as that during which the system 
moves from one essential operating level to another; in this 
case, from one definite value of y to another definite value 
of y. Thus, 


and (18) 


éy(0)= \ 
f'G)=0 


Equation (16) is satisfied only if the integrand is zero, and, 
because ôx 0, the two conditions of equations (16) and 
(17) become 


F, =Å (F) (19a) 


and 

FO)=F jG) 9) (19b) 
Equation (19a) need not hold at t=0 because dy(0)=0. 
The only condition that need hold at ¿=0 is that F;(0) is 
finite, which will be true if y is continuous at f=0. At 
the start of a new transient, 7, F, F,, and F; may be dis- 
continuous, whereas at other points (0<t<4) F; will be 
continuous because of equation (19a). 

Equation (19a) is the differential equation for the y(é) 
that satisfies the original criterion of equation (11). The 
physical answer to the problem is the first integral of equa- 
tion (19a), which satisfies the boundary-condition equation 
(19b). This solution is 


Fy, D=IFiY W (20) 


and whenever y, F;, and so forth are continuous, 


ý |r- Fi) Jo (21) 


by differentiating equation (20) with respect to time. Thus, 


either y=0 (which is true only during stability) or equation 


(19a) is satisfied. 

Thus, equation (20) is the description of that physically 
realizable system the behavior of which will automatically 
and simultaneously satisfy those criterions included in the 
function F during that time interval for which the external 
disturbances are constant and during which the system goes 
from one operating level to any other operating level. A 
stability device must be added to the system; the description 
of such an ideal device is 


y=0 when y=y, (22) 


An additional condition must be met if y is discontinuous 
in the interval 0<é<t,. In this case, F; must be continuous 
during the discontinuity in y; Fp, however, will usually be 
discontinuous when ¥ is discontinuous. This discontinuity 
in F; usually means that some essential external disturbance. 
has entered the system and the point of discontinuity must 
be the start of anew time interval 2 _. 
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APPLICATION TO TURBOJET-ENGINE CONTROLS 


In the usual ease of designing turbojet-engine controls, 
the engine speed that sets the essential operating level of 
the engine and, in the main, other pertinent characteristics 
such as thrust, is to be set or controlled: Limiting conditions 
of the engine are those of. overspeed, overtemperature, 
eompressor surge, and rich burner blow-out. The following 
eriterions on the behavior of this engine are typical: 


at, 
| Fi(N—N, dt for speed control 


« 


ez (N) dt 


for speed overshoot 
for temperature overshoot 


Pir 
Í fs(Dyat and undershoot 
i (23) 


Fi[P—gW)] di for compressor surge 


fd) 
me, 
a 


| P-A WN) at 


fy 
f di 
0) 


Figures 2(a) to 2(e) show the nature of the functions appear- 
ing in expressions (23). The variable involved is essentially 
weighted by the function shown and integrated with respect 
to time. The quantity P—g(N) is the amount the 
compressor-discharge pressure exceeds the safe pressure for 
surge and g(N) (shown in fig. 2(f)) is the compressor-discharge 
pressure for each engine speed at a safe value below surge. 
Rich burner blow-out can be handled in a similar manner 


for blow-out 


for rise time 


(h(N) is shown in fig. 2(g)). The rise time is the total time - 


for the system to move from one essential operating level 
to the other. 

The linearized engine characteristics can be expressed, by 
assuming first-order behavior, as follows: 


T=aN-+aoN (24a) 
P=bN+ceT (24b) 
Thus, the integrals of expressions (23) are of the form 
ti " 
(TAN, Dat (25) 


where f is a continuous function of N and N, and N (barring 
impulsive temperatures and the like) is a continuous function 
of time. f 

Speed control; case A.—If only the error in speed control 
is considered inrportant, the criterion becomes 


f f#(N—N,) di=minimum (26) 


where F=7,(N—N,). Equation (20) becomes 
. A (N—-ND)=0 
and from the nature of f, (fig. 2(@)), (27) 


N=N; 








F N-N;) AN) 
O N-A; No N 
(a) tb} 
A(T) SelP-9(N)] 
Ly Le T g P-giN) 
(c) ( 
fIP-h(W)} g(N) hN) 
PRN) N N 
(e) # (£) (g) 


(a} Speed control. 

(b} Speed overshoot. 

(ec) Temperature orershoot and underskoot. 
(d} Compressor surge. 

(e) Blow-out. 

(f} Compressor line below surge. 

(g) Compressor line below blow-out. 


Figure 2.—Arbitrary weighting functions for various control critecions and pertinent engine 
characteristics involved. 


This result means that, in the absence of other criterions on 
the engine behavior, this speed control should keep speed 
error identically zero, which is physically realizable only in 


the sense of allowing infinite temperatures and the like to — _ 


keep the speed error identically zero. This result, however, 
is inconsistent with the previous development of equation 
(20) in that N is now a discontinuous function of time and the 
time interval of the integral of equation (26) is zero. This 
instance is actually a trivial case of the general problem. 
The result (equations (27)) does indicate that a criterion like 
that of equation (26) must be accompanied by an additional 
criterion (temperature overshoot, for instance) to give a 
physically realizable system. 

Speed control with temperature-limiting criterions; case 
B.—If the error in speed control and the overshoots and 
undershoots in temperature are to be considered the primary 
criterions of control, then from equations (10) and (11), 


if AN-N JHA] dt=minimum (28) 
o 


where F=fi(N—N,)+43(T). Equation (20) then becomes 


NND HMD =M NRT) (29) 


ail 
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and the ideal stability device is such that 
N=0 when N=N, (30) 


Equations (29) and (30) describe the complete control sys- 
tem. 

In figure 2(e), it is convenient to let fy(7)=(T—L,)" for 
Tol: and f,(7)=(T—L;)* for T<L,;. In general, the 
power n should be >1, because, when n<1, T may be infinite 
and of such nature as to make N discontinuous and physically 


É 
unreal eyen though f J(T)dt is finite. In the examples of 


this report, n=2 and f,(N—N,)=(N—N,)*. 

A variety of methods of setting up the control system to 
realize equation (29) exists. When the preceding expressions 
for fı and fs are used, equation (29) can be put in the con- 
venient form 


N) : 
Aa Ney +(L-aNyY=e'oF'N? 
where, for acceleration, when N<N,, 

N>OandL=Z; - 
and, for deceleration, when N>N,, 


N<0 and L=Z, 


(31) 






N 
; ae When N < N, N > O and L*L: 
He~ When N > N,, N < Oond L= L, 
FIGURE 3,—Schematie diagram of control for turbojet engine for case of speed control with 
temperature-limiting criterions (case B). Stability device must be added to this system. 

A schematic diagram of the system is shown in figure 3, 
where equation (31) is considered to give a desired N for any 
value of N, N,, and £. Consistent with equation (31), a 
right-triangle construction is used to give a desired N. The 
actual Ñ, can be made very close to Ñ by using a high-gain 
proportional controller, as shown. Provision must be made 
to change the sign of Ñ and the value of L when N—N, 
changes sign. In addition, the stability condition requires 
a provision for making N=0 whenever speed error is very 
small or zero. 

The use of a high-gain proportional controller, which fol- 
lows from the requirement that N may be discontinuous, 
means that the fuel-flow rate required may be infinite if lags 
exist in the fuel system or in the feedbacks. But, as no 
criterion has been set on fuel-flow rate, this requirement does 
not violate the original criterions. H necessary, however, a 
criterion on fuel-flow rate may then be added to equation 
(28). Even though a criterion on temperature is being satis- 
fied, no direct measurement of temperature is used in figure 3. 
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Actually, the equation for temperature (equation (24a)) is 
used as an indirect indication of temperature. 

The control system of figure 3 has one adjustable param- 
eter X. For any value of A, this system will, for the value 
of integral temperature-overshoot squared obtained, give the 
minimum value of integral speed-error squared. The value 
of \ determines the actual value of the integral temperature- 
overshoot squared. 

The integral temperature-overshoot squared as a func- 
tion of à is shown in figure 4 for the special case where 
aN,=L; that is, acceleration or deceleration to. the speed 
that corresponds to limiting temperature. In this case, the 
system of figure 3 becomes a simple first-order lag system 
and equation (31) becomes 


E(L—aN)=acN 


1 1/2 
E=(1+ a) 


i l Tip T- 

Temperature integral = 3 f TINNA dt 
on 

Maximum temperature = a(N; -NJ j 


Time constant = 2 ae 
ie E Rise time =4 [oy id 
OR 


(32) 


where 








Speed integral = n 





+-K minimum rise time 
for corresponding 
temperature inte- 
gral (from case C 


6 a 4 J 
a ~ clime constant 
a 
if SZ es inimum speed = i 
Dru , integral 


= 


1G 12 14 66 18 20 


r 





F 
























Values of control parameters 






































FIGURE 4.— Various control parameters as functions of Æ for speed control, with temperature- 
Hmiting criterions (case B} when accelerating to limiting temperature. 


In figure 4, the integral speed-error squared, the maximum 
temperature, and the time constant for this transient are 
also shown as functions of à. The curve showing the mini- 
mum rise time for the corresponding temperature integral 
is included for a comparison that will be discussed later. 
The equations for the curves of figure 4 are 





AUTOMATIC CONTROL SYSTEMS SATISFYING CERTAIN GENERAL CRITERIONS ON TRANSIENT BEHAVIOR _ 


; al (GN, aNo nd ae ae 


(pat) dsp 





T L (33) 
+ mar = ee 
ha A 
(time tops ea 
i “og E 


The left sides of these equations have been put in dimension- 
less form. The maximum temperature Tmar occurs at the 
beginning of the transient. The time constant 7 is that for 
the controlled system and is shown compared with the engine- 
time constant s. 

From figure 4, the value of \ is chosen as a compromise 
between the various quantities of equations (33). For H=1 
(A=), the temperature does not overshoot, the speed 
integral is 0.5, and r=. As E increases (à decreasing), 
the temperature integral and the maximum temperature 


aon 











aN 




















FiGCRE 3.—Phase plane showing dyaamics of controlled system for speed control with 
temperature-limiting criterions (case B) where e7A=1, 











increase, whereas the speed integral and the time constant 
decrease. A value of E=,'2 (or a2\=1) appears to be a good 
compromise and is used for the subsequent discussion. 

The behavior of the system of figure 3 can be best seen by 


drawing the trajectories in the phase plane shown in figure. 


5 where acN is plotted against aN for lines of constant 
aN, according to equation (31). Lines of constant temper- 
ature are 45° parallel lines in this plot and the lines of T=L 
are shown. Each trajectory intersects, and is tangent to, 
the line T=L at N=N,. Figure 5 completely describes the 
transient behavior of the system. For any starting point 
anywhere on the plane (for instance, point A), the system 
will automatically move the operating point to that tra- 
jectory corresponding to the N, that exists (point B) and 
then along this trajectory to the point T=L (N=N,, point 
C), and finally the stability condition will enter to move the 
operating point along the solid vertical line to N=0 (point 
D). 

The time sense for these transients is obtained by solving 
the differential equation (equation (31)) for the speed and 
the temperature time responses. The equations for these 
solutions are as follows: 


2u —aN; 


m (GaSe) 14 E— T sinh z (£40) 


EY ( = i)e T sinh z(} +0)+ (34) 


EVEI cosh E (= c)+a — EJ 


T N N 
E 
I Sa 
l ȚsJemperature 
l Ss 
kaana a 
L Se — 
aN; N; 


Engine variable 





£ 


0 ee eT 


a Mo Ng 


FIGURE 6.—Typical transient of controlled system fur speed contro! with temperature- 
limiting criterions (case B) where E='2. 


These transients are shown in the general case in figure 6 for 
a step increase in Ny. Maximum temperature occurs at 
i=0 and the temperature overshoot decreases as NV increases, 
and when N=N,, T=L and T=0. The stability condition 
then causes T to drop to its equilibrium value at this point. 
The time scale shown in figure 6 corresponds to the specific 
relative values of the ordinates shown. 
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Minimum rise time and temperature-limiting criterions; 
case C.—In order to obtain a minimum rise time for a con- 
stant temperature integral, the requirement is that 


f “Emn (35) 


where F=1+A/(7T). Equation (20) becomes 


1+N2(T) = dao Nf (7) (36) 
and the ideal stability device is such that 
N=0 when N=N, (37) 


where N, is the desired stable N. Equations (36) and (37) 
describe the complete control system. Equation (36), and 
therefore the control system proper, does not include N, and, 
in this case, the stability device is independent of the control 
system proper. 

When the same /;(7) is used as was used in the previous 
case, equation (36) can be put in the convenient form 


<+(L—aN)=a%etN 


where, for acceleration, when N<N,, 
N>O and L= L; 

and, for deceleration, when N>N,, 
N<O and L=L, 


(38) 






When N<Ns, N>O ond L-L; 
When N>Ns3, N<O ond L=L, 


FIGURE 7.-Schematie diagram of control for turbojet engine for minimum rise time with 
temperature-limiting criterions (ease C). Stability device must be added to this sytem. 


A schematic diagram of this system satisfying equation’ 


(38) is shown in figure 7. This control system is the same as 


the previous system (fig. 3) except that in figure 7, = 
N—N, : 
replaces —= 
vÀ 
and the stability condition must be imposed outside the 
control system, as N, does not enter into the criterion of 
control (equation (35)). 

The system of figure 7 has one adjustable parameter 4, 
which, as before, sets the precise values of the integrals 
entering into the criterion (equation (35)), as well as all other 
behavior characteristics. 

The temperature integral is shown in figure 8 as a function 
of à for the special case where aN,=Z; that is, for accelera- 


on one leg of the right-triangle construction 
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FIGURE 8.—Various control parameters as function of à for minimum rise time with tempera- 
ture-limiting criterions (case C} when accelerating to limiting temperature. 


tion or deceleration to the speed that corresponds to limiting 
temperature. In addition, the rise time, the maximum 
temperature, and the integral speed-error squared are shown 
for this transient. A curve showing the minimum speed 
integral for the corresponding temperature integral (from case 
B) is also shown and will be discussed later. The equations 
for the curves of figure 8 are as follows: 


=f" ( T—L Ya ¥1+d(L—aN,)? 
=("(-=) a= ie 
ojo \L— aN, (L—a N) À 





1% z +. 
= f di=sinh-' [(L—aN,) yl 


Paz —L J 


ne == 39 
L—aNo Vx (L—aN)) (39) 














(REMY u VTE SNS 
a Jo —aN Q 


2 (L aes aN o) Vn ~ 
l . = 
2x(L—aNo? sinh? afk (L = aN) 


Again à is to be chosen as a compromise between these 


various quantities. In this case, == has the units of 
x 
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temperature and ¥/\(ZL—aN,) is the pertinent dimensionless 
parameter. For a fixed à, the initial speed Na will therefore 
change the time integrals and an increased L—aN, will 
decrease the temperature integral and the maximum tempera- 
ture and increase the rise time and the speed integral. 


adN 










-A =$; L-aNo=3V8 


TÀ = =$, L- -a N= Eves 
Nte 1, L-aNg=ve 











FIGURE 9.~Phase plane showing dynamics of controlled aed for minimum rise time with 
temperature-limiting criterions (case C) where y A (L—aNe} = V3. 


The dynamie behavior of the system of figure 7 7 is shown in 
the phase plane in figure 9. This figure is a plot of aoN 
against a, according to equation (38), for various values of 
à. Lines of constant temperatures are 45° parallel lines and 
the lines T=Z are shown. For any fixed A, only two trajec- 
tories are followed, one for accelerating and one for decelerat- 
ing. From any starting point on the plane, the system will 
automatically move to that trajectory corresponding to ac- 
celeration or deceleration and will move along this trajectory 
until the stability condition enters (at V=N,) to make N=0. 

The dependence of the time mirgral on No may require A 
to vary with No. In figure 9, Ny’s corresponding to each A 
are shown such as to keep NE aN 8. The value of 
this parameter was so taken as to have the temperature 
integral in this case equal to that of the pre ious case for 
purposes of comparison. 








The time sense for these transients is obtained by solving. 


the differential equation (38) for the N—t and T—t tran- 
sients. The equations for these solutions are 


yh (aN—L)=sinh (é +0) 
(40) 
Ee Ge) | 


These transients are shown for the general case in figure 10 
for a step increase in Vy. The temperature will jump to 
some value above L, but in this case, unlike the previous one, 
the temperature continues to increase as the speed increases. 
Maximum temperature occurs at the end of the transient. 
Whenever N=N, (two such conditions are shown), the 


stability condition takes over and T is reduced to its equi- 


librium value. The time scale shown on figure 10 corre- 
sponds to the specific relative values of the ordinates as shown. 
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FısuRe 10.—T ypical transients for controlled system for minimum rise time with temperature- 
limiting criterions (ease C}. 


Comparison of cases B and C.—In case B, a system was 
devised that, for a constant value of the integral temperature- 
overshoot squared, gives a minimum integral speed-error 
squared. In ease C, a system was devised that, for a con- 
stant temperature integral, gives a minimum rise time. For 
case B, the rise time is hot a minimum but can be compared 
with the minimum rise time of case C. For the special 
transient of accelerating to T=L, the system of case B 
reduced to a first-order system, the time constant of which 
is shown in figure 4. Because five time constants are con- 
sidered as the rise time of an exponential, the time constant 
of case B should be compared with one-fifth the minimum 
rise time of case C. In figure 4, for corresponding values of 
the temperature integral, one-fifth the minimum rise time 
from case C is plotted. Figure 4,now shows that the rise 
time for case B is about twice the minimum possible rise time. 


At corresponding values of the temperature integral, 


f” (N—WN,) dé for caseC can be compared with the minimum 
possible value of this quantity from case B. The minimum 


value of 1 (N—N,)? dé is plotted in figure 8. It is seen 
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t 

! (N—N,)? di when 
minimum rise time is obtained is about 120 percent of the 
t 

! (N=N p}? dt. 


that, for the same temperature integral, f 
0 


minimum possible value of f 
JO 


Proof of absolute minimum.—<As previously noted, the 
minimization of equation (11) involves not only the condition 
that the first variation be stationary, which leads to equation 
(20}, but also the condition that the second variation be 
positive. This second condition, however, would prove at 
most a local minimum for the condition of equation (11). 
A special proof of an absolute minimum is shown as follows: 

If equation (11) is written in the form 


f Fo, P EE E E (41) 
Yo yY 


and y is considered as a function of y for the minimum con- 
dition and ġ* is considered as a function of y for any other 
possible case, then the condition for an absolute minimum is 
that 


f ry, ya f "Py pro ug 
vo Y vo Y l 
or 


fiT roi -Ù Fu, i) #>0 (43) 


where dy/y* is a positive differential. If the function F is 
considered in a general form, quadratic in %, as 


FY, YP=GYy) iG W+PE(y) i (44) 
then equation (43) becomes 
Ji@-alrem-@] B20 (45) 


The equation of the control system (equation (20)) becomes, 
for the form of F assumed in equation (44), 


YG, (Y= EY) (46) 


Using this expression for G(y) changes equation (45) to the 
form 


Y 
Í ‘mmMaor—pyco (47) 
Yo y 


For cases B and C, G:=da’o? where A> 0. Thus an absolute 
minimum is proved for these two cases. 

Degree of minimum; case B.—The left side of equation (47) 
can be used as a measure of the deviation from optimum 
conditions when the methods of this report are not used. If 
the two cases (NV for the optimum case and N’* for any other 
case) are compared at the same value of the temperature 
integral, the left side of equation (47) is the difference between 
the integral speed-error squared for any case and the mini- 











mum possible value of this quantity. The ratio of this 
deviation to the minimum value becomes : 


Fractional increase 2E Tt aa dy a 


in speed integral=pe— | ToT 
= r—o 


for the transient of acceleration or deceleration to limiting 
temperature. The coefficient of the integral of equation (48) 
for the value of \ previously chosen (a?\=1) is 2/2. 

Case B can also be considered as giving a minimum value 
of the temperature integral for any definite value of the speed 
integral. If'two cases (Ñ for the optimum and N* for any 
other case) are compared at the same value of the speed 
integral, the left side of equation (45) is proportional to the 
difference between the integral temperature-overshoot 
squared for any case and the minimum possible value of this 





quantity. This deviation can be written 
Fractional increase in 3 EG a) dy (49) 
temperature integral — E —1) Jro — oy* 


for the transient of acceleration or deceleration to limiting 
temperature. The coefficient of the integral of equation 
(49) for the value of \ previously chosen (a?\=1) is 16.4. 
Degree of minimum; case C.—If the two cases N and N* 
are compared at the same value of the temperature integral, 
the left side of equation (47) is the difference between the 


rise time for any case and its minimum possible value. This 
deviation can be written 
Fractional increase A(1—aN) __ z + 
in rise time ~ sinh", (L—aNo) ays z (L—aN o J vo TT 
6 0) 


for the transient of acceleration or deceleration to limiting 
temperature. The coefficient of the integral of equation (50) 
for the value of yA (L—aN ) =2V2 pr eviously chosen is 4.54. 

Case C can also be considered as giving a minimum value 


-ofthe temperature integral for any definite value of rise time. 


If the two cases N and N* are now compared at the same 
value of rise time, the left side of equation (47) is proportional 
to the difference between the integral temperature-overshoot 
squared for any case and the minimum possible value of this 
quantity. This deviation can be written 

= : dy 


vs 
VML—aN, of = oa 


V¥MEL— aNg t 1 pem a 
(51) 





Fractional increase in 
temperature integral — 





for the transient of acceleration or deceleration to limiting 
temperature. The coefficient of the integral of equation (51) 
is 16.4 for the value of 4A (L—aN,) previously chosen. 


GENERAL SUMMARY 


When the criterions on the behavior of a controlled system 
can be expressed in certain general forms, as 
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* 


[Fae 


at 
f ty —y.)’dt 
pn (8) 
(jdt 
eff 


where the time interval is taken as any duration during which 
essential external disturbances are constant and during which 
the system moves from one essential operating point to 
another, the optimum control can be considered as one that 
minimizes one of the integrals of equations (8) while main- 
taining the other integrals at prescribed values. The 
analytical problem, according to the calculus of variations, 
reduces to the following equation: 


fodit | yudit 
ti i 
el Flodin f dt=minimum (9) 


For general first-order systems, where z==z(y,y), equation 
(9) reduces to 


t 
[ ‘Bay dt =n (11) 


where F=f(y)--Mi(y—ys)?+rAfo(z) +3. The equation nec- 
essary for the control system to satisfy equation (11) and 
all the boundary conditions becomes 


FY, N= 9Fi yy (20) 


This equation should be followed by the control system 
proper. In addition, a stability device must be added to the 
system, the idealized characteristics of which would make 
y=0 when y=y,. 

The arbitrary multipliers à are then found by evaluating 
the integral criterions involved in F. The transient behavior 
of the system derived is.found by solving the differential 
equation (20). The degree of the minimum or the amount 
suffered when any other control system is used was evaluated 
for the special cases considered. A summary of these devel- 
opments follows for a special form of the F function, quad- 
ratic In 47, where 


Fly, n=@y)+yGity) +y@s(y) (44) 


Control-system equation.—For F in the form of equation 
(44), the control-system equation (20) becomes 





«> Gy) 
Y — Gy) (32) 


The function G,(y) does not affect the control system and 


272483—54——_15 








the control-system equation gives an explicit expression for 
y as a function of y. 

Evaluation of integrals—If an integral of a function 
H(y,y) is any one of the criterions to be considered, it can be 
evaluated as follows: 


t E b 2, Gy) | (GM) 4 eo 
f Heyindim fH] a4 VEI VT © 


Thus the integrals can be evaluated without solving the 
differential equation (52). 

Transient behavior.—The differential equation (52) can 
be easily solved as follows: 


Giz) _ ~ 
ENS Ga) dx=t—ty (54) 


Degree of minimum.—Equation (47) was derived to prove 
the absolute minimum for the special cases considered and 
can be used to evaluate the degree of the minimum found. 
Thus, 











it TE ty 3 a Seat iia dy 
f, Aly, y*)dt f Hy, ydt= | COB AIA) FO 
(55) 
where i A(y,y)di is to be a minimum or maximum. 


If GYO for y»<y<yy,, then an absolute minimum is 
obtained; if Galy) <0 for yo<y<y;, an absolute maximum is 
obtained. Equation (55) also indicates that the degree of 
minimum (or maximum) varies with the magnitude of G,(y). 
The y* in the denominator indicates that for any deviation 
ly*—y| in ġ, it is better to err on the side of a larger |y*1. 


SUMMARY OF RESULTS 


A rational analytic method: for the design of automatic 
control systems was derived. Criterions on the behavior of 
the controlled system were developed in the form of certain 
time integrals. When any of these arbitrary but physically 
realizable criterions were used as a starting point, those 
equations that the control system must follow were derived. 
The criterions developed required the minimization of certain 
time integrals using the calculus of variations. The method 
gave not only a description of the behavior of the controlled 
system but also those physically realizable equations that the 
eontrol system can follow in order to realize this behavior. 
General results were shown for systems of second order and 
of any number of degrees of freedom. 


Lewis ELIGET PROPULSION LABORATORY 
NATIONAL ÅDVISORY COMMITTEE for AERONAUTICS 
CLEVELAND, Onto, October 11, 1950 
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APPENDIX ' 


SECOND-ORDER SYSTEMS OF SEVERAL DEGREES OF FREEDOM 


It is beyond the scope of this report to detail the cases of 
higher-order systems and those having several independent 
variables. The general equations for these cases will be 
developed and it is expected that subsequent reports will 
cover their applications in detail. 

For the case of a second-order system with two degrees of 
freedom, equation (11) becomes 


ie F(y,3,4, 2,4, 3) dt=minimum (A1) 


The 


where y and z are independent functions of time. 
condition to satisfy equation (A1) is 


d ty +edt; . n , 
del,  FOteyýtesjütesje+esz, 
2+662,2+ €62)di=0 


(A2) 
at 
e=0 


The time duration of the integral of equation (A1) is that 
beginning at a definite time (=0) but not ending at any 
definite time but rather along some curves: y=fi(t), y=h(t), 
z=g.(t), and 2=g,(t) (fig. 11). The functions ôy and ôz are 
arbitrary and independent functions of time. 
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Ficur& 11,—Hlustration of curves for minimization of f Fiy, ýt, 2) dt. 
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Performing the operation indicated by equation (A2) 
gives 


f "P 5y-+F)8)-+ Fjot Foz -+ F, t Fiddi +F E= 0 


(A3) 
After integration by parts, 


f [r-i roti E ora [RF eat 
Ta P) | 2 ate (Fi ay-+Frad—F, Fiy) | rE 


sg f 
(Mise Fi; (Fase) | =0 (A4) 


As before, the integrands of the integrals and the boundary- 
condition terms must vanish separately. From the geometry 
at the end condition (f=¢,, fig. 11), 


òy =fr E) iE), t: 
y= [f 6) GE) bh: 


l (45) 
ôzt)= fg E)E t), oe, 
òi) = [g hh) — 2] t 
The three conditions from equation (A4) then become 
d d? 
F- Ci+gp (Fy) =0 
d d 
F- d+ 7p (F)=0 
FoF, +94 FoF iF 22(F)—aF | d 
Yi TY g UYE a ave ae a leek, 
a l , 
O| rig E | + HG P+ (A6) 


alt) | P-gp E | tee Ft 
wO | Fig En |, O Fo 
8200) | Fi— i; E» | +220) F0)=0 | 


The first two of equations (A6) are the differential equa- 
tions that satisfy the original criterions of equation (A1). 
The physical solution to the problem is the pair of solutions 
of these equations that satisfies the boundary-condition 
equation of equations (A6). If the first of equations (A6) 
is multiplied by y and the second, by 2 and the equations are 
added, an exact derivative is formed, the integral of which 
is 


FiFi +9 5 0) oF F428 E EFFO (AN) 
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A study of the boundary-condition equation of equations 
(A6) then shows that the physically reasonable boundary 
conditions should be as follows: 


If Fy #0 or Fy AK, then sy(0)=0 and f,/¢)=0 
If Fz 0, then 67(0)=0 and f,’(f,)=0 
If Ff; 0 and F;#K, then 62(0)=—0 and g/'¢,)—0 
If Fs0, then 62(0)=0 and g,’(,)=0 


(A8) 


Then in equation (A7), C=0 and the final solution to the 
problem of equation (A1) for the boundary conditions of 
equations (A8) is 


Fog Fi +9 5 Eng Fy —2F 42 & F)—2F <0 


and either 


F, -4 Pit Fy=0 + (Ad) 


or 


F= F+ F)=0 


The meaning of the boundary-condition equations (A8) 
is to define the original criterions for that duration during 
whieh the system moves from one essential operating level 
to another. Thus, if all conditions of equations (A8) must 
hold, the system goes from one definite y, 7, z, and Zz to any 
other definite y, 4, z, and 2. The first differential equation 
of equations (A9) would be of third order and the second or 
third equations of (A9) would be of fourth order. If equa- 
tions (AQ) are integrated to obtain a pair of second-order 
differential equations having three constants of integration, 
the choice of these three constants can then determine a 
desired end point, that is, the values of ý, z, and Z at some 
final y. 

The physical solution to the problem, then, is the pair of 
second-order differential equations that are solutions of 
equations (AQ) with the constants of integration evaluated 
so that the system goes through some desired end point. 
This end point must be such as to allow the possibility of 
stability. Such an end point may be written as follows: 








y=0, z=2,, and 2=0 when y=y, (A10) 


which gives three conditions for the evaluation of the three 
constants of integration. A stability device must still be 
added to the system so that, at the point when equations 
(A10) hold, 


ğ=0 
and l 


=0 


(A11) 


Mir 


Equations (A9) are symmetric in the variables y and z, 
which indicates the nature of the extension for more inde- 
pendent variables. Thus, if a third independent function 
w exists, the original criterion would be written 


(A12) 


| F (y, ý, 9, 2,2, 2, w, 8, W) dt=minimum 


This condition is satisfied under boundary conditions similar 
to equations (A8), where two additional conditions are added 
on the variable w,and the following equations deseribe the 
controlled system: _ 


‘ .d 2 TOPET : 
F—ýFi+ý q Pift iti g EdF 
ES, a ee ee 
wher wb a (Fy)—vF g= 
and any two of the following three equations: 


d a (A13) 
Mya Pit ap (F) =0 


d d? Ap 
Fe (F; tJe (Fj=0 
d epa @ opa 
Pua Pat 7p F= 


The three equations of equations (A13) can then be integrated 
to give three second-order differential equations where the 
five constants of integration are evaluated so that the system 
goes through the desired value of ġ, z, 2, w, w for some final 
value of y. 
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